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ABSTRACT 

The distribution of galaxies displays anisotropy on different scales and it is often referred as galaxy 
alignment. To understand the origin of galaxy alignments on small scales, one must investigate how 
galaxies were accreted in the early universe and quantify their primordial anisotropic at the time of 
accretion. In this paper we use N-body simulations to investigate the accretion of subhaloes, focusing 
on their alignment with halo shape and the orientation of mass distribution on large scale, defined 
using the hessian matrix of the density field. The large/small (el/e3) eigenvalues of the hessian matrix 
define the fast/slow collapse direction of matter on large scale. We find that: 1) the halo major axis 
is well aligned with the e3 (slow collapse) direction, and it is stronger for massive haloes; 2) subhaloes 
are predominately accreted along the major axis of the host halo, and the alignment increases with 
the host halo mass. Most importantly, this alignment is universal; 3) accretion of subhaloes with 
respect to the e3 direction is not universal. In massive haloes, subhaloes are accreted along the e3 
(even stronger than the alignment with the halo major axis), but in low-mass haloes subhaloes are 
accreted perpendicular to e3. The transit mass is lower at high redshift. The last result well explains 
the puzzled correlation (both in recent observations and simulations ) that massive galaxies/haloes 
have their spin perpendicular to the filament, and the spin of low-mass galaxies/haloes is slightly 
aligned with the filament, under the assumption that the orbital angular momentum of subhaloes is 
converted to halo spin. 

Subject headings: cosmology: theory - dark matter - large-scales structure of Universe - galaxies: 
halos - methods: statistical 

1. INTRODUCTION 

In the cold dark matter universe (CDM), structures 
emerge from the initial seed of perturbations via gravita¬ 
tional instability. On large scales, the mass distribution 
is characterized as cosmic web (filaments, sheets, clus¬ 
ters, voids), which can be fairly described by a linear the¬ 
ory and the Zel’dovich approximation (Zel’dovich 1970). 

On small scales, dark matter haloes form by collapse of 
matter at the intersection of filaments, and they gradu¬ 
ally merge and grow in a hierarchical and non-linear way. 

The detail of the dark matter halo properties and mass 
distribution on large scales can be accurately studied by 
N-body simulations (e.g. Springel et al. 2006). Under 
such a formation scenario, it is naturally expected that 
the properties of haloes are correlated with the large- 
scale structure (LSS). For example, the halo shape is well 
aligned with the cosmic filament (e.g., van Haarlen & van 
de Weygaert 1993; Hahn et al. 2007a; Aragon-Calvo et 
al. 2007a; Zhang et al. 2009; for a summary, see Forero- 
Romero et al. 2014 ), the halo spin is originated from the 
tidal torque force by the mass distribution on large scales 
(White 1984), and its direction is well correlated with the 
cosmic web (e.g., van de Weygaert & Bertschinger 1996; 

Bond et al. 1996; Lee 2004; Bailin & Steinmetz 2005; 

Hahn et al. 2007a; Zhang et al. 2009; Libeskind et al. 

2012). Most importantly, the alignment between halo 
spin and the LSS has a dependence on halo mass (e.g., 

Aragon-Calvo et al. 2007a; Hahn et al. 2007b; Codis 
et al. 2012; Trowland et al. 2013; Dubois et al. 2014; 
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Zhang et al. 2015). 

These theoretically predicted correlations of the halo 
properties with LSS are manifested by the distribution of 
galaxies inside dark matter haloes, as galaxies are ’lumi¬ 
nous’ tracers of merged haloes formed at early times. For 
example, satellite galaxies are found to be aligned with 
the major axis of central galaxies, with strong depen¬ 
dence on galaxy colors and halo mass (Brainerd 2005; 
Yang et al. 2006). Local dwarf satellites in the Milky 
Way and M31 are distributed anisotropically and most 
are in a large thin disk (Kroupa et al. 2005; Metz et 
al. 2007; Ibata et al. 2013). The shapes of luminous 
red galaxies are correlated with each other up to very 
large scales (e.g., Faltenbachler et al. 2007; Okumura 
et al. 2009; Li et al. 2013). Galaxy spin is also found 
to correlate with LSS although the measurements of the 
galaxy spin are not straightforward (e.g., Trujillo et al. 
2006; Jones et al. 2010). More interestingly, it has been 
recently found that the correlation of the galaxy spin 
with the LSS depends on galaxy type/mass, such that 
the spin of early type galaxies is perpendicular to the 
filament around, but late-type galaxies have their spin 
slightly parallel to the filament (Tempel & Libeskind 
2013; Zhang et al. 2015), in amazing good agreement 
with the theoretical expectations (e.g., Aragon-Calvo et 
al. 2007a). 

The observed anisotropic distribution of galaxies on 
different scales raises the question of its origin. Some 
studies (e.g., Benson et al. 2005; Wang et al. 2005; 
Libeskind et al. 2005) have shown that satellite galax¬ 
ies are accreted anisotropically during the formation of 


2 


dark matter haloes, and this primordial anisotropy is the 
origin of the current anisotropic distribution of satellite 
galaxies. However, the observed distribution of galax¬ 
ies within DM haloes involves complicated baryonic pro¬ 
cesses (e.g., Kang et al. 2005; Vogelsberger et al. 2014) 
and galaxies have experienced the non-linear dynamical 
evolution inside dark matter haloes. Along with the fact 
that the signal of the galaxy anisotropic distribution is 
much weaker, it is therefore difficult to infer how much of 
the currently observed anisotropic distribution of galax¬ 
ies is from the primordial anisotropy set on large scales 
before galaxies were accreted, and how much of it is from 
the non-linear evolution inside the dark matter halo. For 
example, the anisotropic distribution of satellites respect 
to central galaxies can be purely ascribed to the non¬ 
linear evolution or the non-spherical nature of dark mat¬ 
ter haloes (e.g.,Jing & Suto 2002; Kang et al. 2007; Wang 
et al., 2014; Dong et al. 2014; Debattista et al. 2015). 
For a recent review on galaxy alignments and its rela¬ 
tion with the halo shape or the LSS, see Joachimi et al. 
(2015). The flip of the halo spin-LSS correlation (e.g. 
Aragon-Calvo et al. 2007a) also attracts great atten¬ 
tion recently. Using hydrodynamical simulations, it was 
found that the flip is related to the halo merger history 
and cold gas accretion during the formation of elliptical 
and spiral galaxies (e.g., Codis et al. 2012; Welker et al. 
2014; Codis et al. 2015). 

To understand the observed correlation of the galaxy 
distribution with respect to the halo shape or the LSS 
in detail, one needs to trace galaxies back to the time of 
their accretion, to quantify the degree of the primordial 
correlation with DM haloes or the LSS at that time. This 
is often done by tracing subhaloes back to early times us¬ 
ing N-body simulations. In this paper we study the spa¬ 
tial alignment of subhaloes with respect to the host halo 
shape and the LSS, and in a future paper we will focus 
on the velocity/spin correlation with LSS at the time 
of accretion. There were only a few studies using this 
kind of approach. Benson et al. (2005) found that sub¬ 
haloes are preferentially accreted in the planes defined by 
the major and middle axes of the host haloes. Wang et 
al. (2005) also found that subhaloes are accreted along 
halo major axis and massive subhaloes show a stronger 
trend. However, they both ignored the correlation with 
the LSS. Recently, Libeskind et al. (2014) and Shi et al. 
(2015) find that subhaloes are predominantly accreted 
along the filament, but their classification of cosmic web 
is based on velocity shear or tidal field. In our paper we 
classify the cosmic web by means of the density field in 
order to classify the LSS environment, and our analysis 
goes to very low-mass haloes. In particular, we will show 
that our results are helpful to explain the puzzled non¬ 
universality in the correlation between halo/galaxy spin 
with the LSS. 

The paper is organized as follows. In Section.2 we in¬ 
troduce the simulations and the methods to quantify the 
halo shape and the LSS. In section.3 we show how the 
halo shape is correlated with the LSS and the anisotropic 
distribution of accreted subhaloes with respect to the 
halo shape and the LSS. Then, we present their depen¬ 
dence on halo mass and redshift. In section.4 we summa¬ 
rize our results and briefly discuss how our results help 
to explain the correlation between the halo spin and the 
LSS. 


2. N-BODY SIMULATION AND LSS CLASSIFICATION 

In this work we use two N-body simulations with differ¬ 
ent box sizes. Both simulations are run with the Gadget- 
2 code (Springel 2005) using the same number of par¬ 
ticles of 1024^ and the same cosmological parameters 
from the WMAP7 data (Komatsu et al. 2011, namely: 
Da = 0.73, Dm = 0.27, h = 0.7 and erg = 0.81. The 
box sizes are 200Mpc/h and 65Mpc/h, and the particle 
masses in these two simulations are 5.5 x 10^Mq/H, and 
1.8 X IO^Mq/Zi, respectively. With the larger simula¬ 
tion box we get more well resolved massive haloes, while 
the smaller box one allows us to resolve haloes down to 
low masses, about 5 x 10^Mq/H. We note that in the 
following analysis the results for haloes with mass be¬ 
low IO^^Mq are from the simulation with small box size 
{65Mpc/h), and for massive haloes we use the results 
from the 200Mpc/h simulation. 

From the simulation output we firstly identify dark 
matter haloes using the standard friend-of-friend (FOF) 
algorithm, and each FOF halo should contain at least 20 
particles. Then, subhaloes in each FOF halo are found 
using the SUBFIND algorithm (Springel et al. 2001) 
and their merger trees are also constructed. For details 
the readers can refer to Kang et al. (2005). We trace 
each subhalo back to the time when it was last identified 
as an independent FOF halo (labeled as h). In a very 
short time period (the next snapshot) the small halo h 
will merge with a larger halo (host halo, labeled as H). 
By means of positions and velocities of the haloes h and 
H, we can obtain the time and position when the halo 
h crosses the virial radius of halo H for the first time. 
Note that in the following, when we refer to the subhalo 
distribution at accretion, we are actually referring to ei¬ 
ther the distribution of the halo h with respect to the 
host halo H, or with the large scale around the halo H 
at the time when halo h crosses the virial radius of H for 
the first time. 

For each FOF halo its virial mass is the enclosed mass 
within the virial radius inside which the mean density is 
Ac{z) times the critical density of the Universe (Bryan 
& Norman 1998). To measure the shape of each halo, we 
follow the traditional method to calculate the normalized 
inertia tensor defined as (Bailin & Steinmetz 2005), 

= ( 1 ) 

where a;q„ is the distance component of the particle n to 
the halo center, and the summation is over all particles 
inside the virial radius. The three eigenvectors of lij de¬ 
fine the orientation of the three axes of a halo, and the 
direction of the largest eigenvalue defines the major axis. 
Note that the minimum number of particles in each host 
halo H is taken as 500, so the halo shape is more accu¬ 
rately determined (Jing & Suto 2002). We also tested 
the halo shapes (the three eigenvectoes) using all parti¬ 
cles of the FOF halo, and found that they agree quite 
good with those from those using the particles inside the 
virial radius. 

To define the LSS environment around each halo, we 
follow the most used method (e.g., Hahn et al. 2007a; 
Arag’on-Calvo et al. 2007b; Zhang et al. 2009,) by cal¬ 
culating the Hessian matrix at the position of each halo. 
The Hessian matrix of the smoothed mass density field 
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Fig. 1.— Slices from the simulation. Left: the 2D mass distribution around a halo at z = 1, the red circle denotes its virial radius, and 
blue circles denote those subhaloes to be merged with the host halo in a short time (some subhaloes are located inside the red circle due 
to projection effect). Yellow and cyan lines show the e3 direction of the density field on large scales, and the major axis of the halo. Right: 
the slice at 2 = 0, the halo is the descendent of that shown in the left panel, and blue circles show the z = 0 distributions of those accreted 
subhaloes labelled in the left panel. 


is defined as, 


, ^ d'^Psjx) 

dxidxj 


( 2 ) 


where Psix) is the smoothed density with smooth length 
Xs- There is lack of systematic census on which smooth 
length is the best to characterize the mass distribution 
on large scales (for partial discussions, see Hahn et al. 
2007b; Forero-Romero et al. 2014). In most studies, a 
constant smooth length (0.5 ~ 2Mpc/h) is used at z = 0 
(e.g., Hahn et al. 2007a; Zhang et al. 2009; Codis et al. 
2012; Trowland et al. 2013). However, it is found that 
an evolving smooth length is more physical reasonable 
and can better describe the evolution of the LSS envi¬ 
ronment (e.g., Hahn et al. 2007b; Libeskind et al. 2014). 
Following Hahn et al. (2007b), we adopt a smooth length 
roughly scaled as the virial radius of the typical mass for 
halo collapse at different redshfit (often written as M*). 
To be able to compare with other studies at z = 0, we 
adopt Ts ^ 2/{l-\-z)Mpc/h. Actually, we also found that 
our results are not significantly affected by using other 
constant smooth length (1, 2, 5 Mpc/h) at all redshifts. 
We will show some comparison in Fig. [ 2 ] 

The eigenvalues (Ai > A 2 > A 3 ) of the Hessian ma¬ 
trix define the LSS environment around a halo, and it 
is classified as cluster, filament, sheet and void depend¬ 
ing on the number of positive eigenvalues. For example, 
there is no positive eigenvalue for a cluster, one positive 
eigenvalue for a filament, and two positives for a sheet. 
This classification of halo environment well mimics the 
description of Zel’dovich theory (for a review, see Cautun 
et al. 2014). Basically, the eigenvector ei (corresponding 
to the largest eigenvalue Ai) defines the fast collapse di¬ 


rection, and 63 corresponds to the slowest collapse direc¬ 
tion. For example, for a halo in a filament environment, 
63 gives the direction where the matter has not collapsed 
on large scales, and future mass accretion will mostly 
happen along 63 . As the matter distribution along 63 
has not strongly collapsed, the halo is less compressed 
along this direction, and this is why the halo major axis 
will then have a tendency to align with 63 , which has 
been shown in many studies (e.g., van Harrlem & van de 
Weygaert 1993). 

Fig. [ 1 ] illustrates the space configuration described 
above. The left panel shows the mass distribution around 
a selected host halo at z = 1, with virial radius labelled 
as the red circle. The blue circles denote the subhaloes 
which will merge with the host halo in a short time scales. 
The yellow line shows the 63 direction, which is seen to be 
well aligned with the mass distribution on large scales. 
Cyan line is the direction along the major axis of the 
host halo. The position angles of subhaloes with the 
cyan/yellow line are 9mA and respectively. The right 
panel shows the distribution of those accreted subhalo at 
z = 0, their anisotropic distribution with respect to the 
host major axis can also be seen in the plot. 

With the above descriptions, we are able to derive the 
spatial distribution of subhaloes with respect to the host 
halo and its LSS environment around. The angular po¬ 
sition of the halo h at the time of accretion with respect 
to the major axis of the halo H is labeled as and 
the angle respect to the slowest collapse direction ( 03 ) 
centered on H is labeled as For a random distri¬ 
bution, the expected value of < |cos(d)| > is 0.5, and 
if < I cos (0)1 > is larger than 0.5, we refer to it as an 
alignment with the halo major axis or the LSS. 
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Fig. 2.— The mean angle between the halo major axis and the slowest collapse direction on large scales, eg, calculated from the Hessian 
matrix using different smooth length. The left panels show the dependence of the alignment on the smooth length at given redshift, and 
the right panels show the redshift evolution of the alignment for a given smooth length. 


3. RESULTS 

3.1. Accretion along halo major axis and the LSS 

Observations of galaxy distributions have found that 
satellite galaxies are aligned with the major axes of their 
host galaxies, and this alignment is dependent on the 
mass/color of the hosts (e.g., Yang et al. 2006). Due 
to the fact that the observed signal is weak compared to 
that from N-body simulations (e.g., Kang et al. 2007), 
satellite galaxies alignments can be explained if the cen¬ 
tral galaxy roughly follows the shape of the dark matter 
halo and the tracing becomes stronger with halo mass 
(e.g., Dong et al. 2014). Such an explanation implies 
that the observed alignment can be purely ascribed to 
the non-linear evolution inside DM haloes. 

In the first part of this section we study if such an 
alignment between subhaloes positions and halo major 
axis (or eg of the LSS) can be seen at the time of accre¬ 
tion, and investigate which kind of alignment is stronger 
{Oma or 6 * 63 ). Before discussing our results, we first check 
how the shape of dark matter haloes is correlated with 
the tidal field around them. In Fig. we plot the av¬ 
erage angle between the major axis of the halo and 63 , 
the slowest collapse direction of the mass on LSS. In the 
left panel, we show the dependence of the alignment on 
smooth length (1, 2, 5Mpc/h), and in the right panel we 
show the redshift evolution for a fixed smooth length. 


Note that here the halo mass is scaled by the character¬ 
istic mass (M») of halo collapsing at given different. By 
scaling the halo mass with M* we are actually studying 
the peak of the density field (with given height) with its 
LSS around, which has a more physical meaning (Hahn 
et al. 2007b; Trowland et al. 2013;) 

Fig.dlshows that the major axis of the DM halo is well 
aligned with its LSS ( 03 ), and the alignment increases 
with halo mass. The left panel shows that at given 
redshift, the alignment is stronger for a smaller smooth 
length (IMpc/h), indicating that the mass distribution 
at large scale is better described by this smooth length. 
The right panel shows that at a given smooth length, the 
alignment is evolving with redshift, and being stronger 
at low redshift due to the non-linear evolution. The uni¬ 
versal alignment between halo major axis with the orien¬ 
tation of the density field {eg) agrees well with the pub¬ 
lished results so far (e.g., Forero-Romero et al. 2014, and 
references therein). The results clearly show that when 
a halo collapses its internal mass will be strongly com¬ 
pressed along the fast collapse direction where the mass 
over-density is highest, so the particles will more likely 
move along the least compressed direction eg. Such a 
scenario well captures the spirit of Zel’dovich theory and 
is a manifestation of the correlation between halo shape 
(density peak) and the orientation of the density field on 
large scale (van de Weygaert & Bertschinger 1996; Bond 
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Fig. 3.— The alignment of subhaloes with the halo major axis and the es direction of the LSS at their time of infall (see text for detail). 
Different panels are for the final host haloes at 2 = 0 with different mass. The alignment with the halo shape (black solid lines) is universal 
and increases with halo mass. The alignment with 63 (red dashed lines) is stronger in massive haloes, but it is reversed for low-mass haloes 
where it is perpendicular to 63 . 


et al. 1996; Lee & Pen 2000; Porciani et al. 2002; Rossi 
2013; Libeskind et al. 2014). 

Now we investigate the primordial anisotropy of sub¬ 
haloes at the time of accretion. As found in observa¬ 
tions, the alignment of satellites is stronger in massive 
host haloes. For this reason, we investigate if the pri¬ 
mordial anisotropy is also dependent on the mass of host 
haloes. We select host haloes at z = 0 with a wide range 
of mass from IO^^Mq to IO^^Mq, and trace their sub¬ 
haloes back to the accretion times. For each accretion 
event, we can get the angle between the infalling po¬ 
sition of a subhalo (halo h) and the major axis of the 
host halo, cos{0ma), and the angle with the 63 direc¬ 
tion, cos{9f,^), respectively. The distributions of them 
are shown in Fig. [31 and different panels are for host 
haloes with different virial mass at z = 0. The black 
solid lines are distributions of \cos(9ma)\ and red dashed 
lines are for [cos( 6 * 63 )!. 

Some interesting results can be seen from Fig. [31 
Firstly, the infalling subhaloes have a tendency to align 
with halo major axes, and the alignment increases with 
host halo mass. However, the alignment with 63 is differ¬ 
ent and, more importantly, it is not universal. In massive 
haloes (lower right panel), the alignment is even stronger 
than the halo major axis, but it is anti-aligned in low- 


mass host haloes (upper left panels), where subhaloes 
are accreted perpendicularly to 63 . So, the accretion in 
low-mass haloes is preferentially along the fast collapse 
direction. This has important implication in explaining 
why the spin of low-mass haloes is aligned with the fila¬ 
ment. We will show the results in detail in Section. 13.21 
The Zel’dovich theory predicts that structure forma¬ 
tion proceeds subsequently along the three eigenvectors 
of LSS, i.e, it firstly collapse along ei to form a sheet, 
and then the second collapse happens along 62 to form 
a filament. A halo will form if the matter along the fil¬ 
ament ( 63 ) continues to collapse. In this secnario, the 
accretion of a growing halo is expected to along the re¬ 
cent collapsed direction ( 63 ). Although the Zel’dovich 
approximation is only valid on linear scale in Lagrangian 
space, it remains valid to describe the feature in the Eu- 
lerian non-linear regime, as we see in the lower panels. 
Only on very small scales the accretion departs from the 
predictions of the Zel’dovich theory (upper left panels). 
On small scales, the mass accretion is more randomly, 
and being slightly aligned with the halo major axis, im¬ 
plying that mass accretion is more likely determined by 
the local halo potential. Departure of anisotropic accre¬ 
tion from the Zel’dovich prediction has been realized for a 
long time (e.g., Icke 1973; Silk & White 1978; Eisenstein 
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Fig. 4.— The current alignment of subhaloes with the halo major axis at z=0 (red dashed lines), and the alignment at accretion (black 
solid lines). It is found that the current alignment with the halo major axis is higher than that in the past, indicating a non-linear evolution 
effect which is more obvious in intermediate halo masses. 


& Loeb 1996). Bond & Myers (1996) provided a compre¬ 
hensive description about the nonlinear anisotropic col¬ 
lapse of ellipsoids, which matches better the simulation 
results than the Zel’dovich theory. 

Our results about subhalo accretion with ea marginally 
agree with the recent results of Libeskind et al. (2014) 
and Shi et al. (2015), who have also studied the accre¬ 
tion of subhaloes with respect to the LSS. Libeskind et 
al. (2014) found that subhaloes are predominantly ac¬ 
creted along the 63 direction, and this effect increases 
with both the subhalo and host halo masses. As we will 
see in the next section, they obtained a positive align¬ 
ment because their bins of host halo mass is wider, thus 
the anti-alignment from low-mass haloes is immersed in 
the positive alignment signal from high-mass haloes. Shi 
et al. (2015) use a method similar to ours to define the 
large scale structure, and they also found that subhaloes 
are accreted along the 63 direction. However, their anal¬ 
ysis is only for massive haloes (> 10^'^Mq) where their 
results are consistent with ours. It is not clear what are 
the results for subhalo accretion in the low-mass haloes 
in Shi et al. (2015). 

In Fig. 0] we compare the current alignment of sub¬ 
haloes (using subhaloes’ positions and the host halo ma¬ 
jor axes at z = 0) with that at infall. The black solid 
lines are the same as in Fig. [3l and the red dashed ones 


show the current alignment in different host haloes. It 
is seen that the z = 0 alignment is stronger in massive 
haloes, which is consistent with what the data suggests. 
In general, the current alignment of subhaloes with their 
hosts is larger than the alignment at the time of accre¬ 
tion. A stronger evolution is more obvious in haloes with 
mass around IO^^Mq. However, observed galaxies often 
live in haloes with mass larger than IO^^Mq (Yang et al. 
2006), and the observed alignment of galaxies is much 
weaker, even compared with the primordial alignment 
(black solid lines). So, it is still difficult to quantify how 
much of the observed signal is from the primordial align¬ 
ment and from the evolution. 

In Fig. [3] the distribution of accretion angle is a pure 
number counting, and not weighted by the mass of sub¬ 
haloes. In Fig. Owe plot the cumulative fraction of mass 
accreted along the 63 direction (left panel) and along the 
major axes of haloes (right panel). Note that in our anal¬ 
ysis the minimum halo mass counts only 20 particles, so 
we neglect the smooth mass accretion and consider only 
the resolved mass in subhaloes. It is seen that in mas¬ 
sive haloes (> 10^'^Mq), the net mass accretion is along 
the major axis of the host or along 63 , and the latter is 
slightly stronger. For example, more than 50% of the re¬ 
solved mass is accreted with 9 < 30°. For low-mass host 
haloes, the mass accretion is close to be isotropic. 
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Fig. 5 . — The cumulative fraction of accreted mass in subhaloes with a given alignment angle. Left panel is for alignment with 63, and 
right for alignment with halo major axis. A higher fraction of mass is accreted along 63 in massive haloes, and it becomes almost isotropic 
in low-mass haloes. 



Fig. 6 .— The alignment of accreted subhaloes with the host halo major axis (left panel) and the 63 direction of the LSS. Different 
from previous figures, here no constraints are posed on the final host halo mass, and the accretion angles are binned depending on the 
accretion events with different host halo mass and redshift. As we can see, the accretion is always along the halo major axis (left panel). 
The alignment with 63 is not universal and it is perpendicular to £3 for low-mass haloes, and the flip mass is lower at high redshift (right 
panel). 


3 . 2 . Dependence on accretion mass and redshift 

Above we have shown the overall accretion of subhaloes 
along the halo major axis and 63 of the LSS. However, 
in that analysis we fix the host halo mass at 2 ; = 0 and 
study the accretion of their progenitors at high redshift. 
Doing that is to understand the primordial alignment 
of accreted subhaloes and to compare with their current 
alignment in the z = 0 hosts. In this case, the main 
progenitors and the accreted subhalo mass are all evolved 
with redshift. In this part, we do not fix the z = 0 
host, but study the alignment for fixed host halo mass at 
different redshift. In this way, we can quantify the mass 
dependence and being able to compare with the recent 
study by Libeskind et al. (2014). 

In Fig. | 6 ] we show the alignment along the host halo 
major axis (left panel) and the 63 direction (right panel). 
The left panel shows that the alignment with halo major 
axis increases with host halo mass, and it is stronger 
for haloes at high redshift. Note that here we do not 
scale the halo mass by the characteristic mass of halo 
collapsing (M*), as we did in Fig. [H So for given halo 
mass, the density contrast is higher at high redshift. Note 


that for massive haloes the alignment from the simulation 
of 200Mpc/h box is systematically lower than that from 
the small box simulation, and this is due to some low- 
mass subhaloes that have not been resolved in the large- 
box simulation. Basically, our results agree with previous 
studies (e.g., Benson 2005; Wang et al. 2005; Shi et al. 
2015) which state that the alignment with the halo major 
axis is stronger for massive haloes. It is also seen that 
the alignment is universal, i.e., that subhaloes are always 
accreted along the major axis although it is close to be 
isotropic for low-mass haloes. 

The results in the right panel are more interesting. 
First, the accretion with 63 of the LSS is stronger than 
the alignment with the halo major axis shown in the left 
panel, and it also increases with host halo mass and red¬ 
shift. Moreover, such dependences on halo mass and ac¬ 
cretion time are much stronger than the alignment with 
the halo major axis. Second, it is seen that the alignment 
with the 63 is not universal. For low-mass haloes, th ac¬ 
cretion is actually perpendicular to 63 , and the transit 
mass of the alignment is around 3 x at z = 0, 

and it decreases to 5 x IO^Mq/H at z = 2. 





















Libeskind et al. (2014) also studied the alignment 
of subhaloes with the 63 axis of the LSS around host 
haloes. However, their approaches are different from 
ours. Firstly, they use the velocity shear to define LSS, 
whereas we use the density field (see Forero-Romero et al. 
(2014) for discussion and comparison between the classi- 
hcation of the LSS based on density and velocity field). 
Secondly, the alignment of subhaloes in their analysis 
are measured when subhaloes cross a few virial radius 
(4 ~ 16i?„ir) of the host halo. Although the approaches 
are quite different, we get similar dependences on host 
halo mass and redshift compared with their results. The 
main difference is that they found that the accretion is 
always along 63 , while we find that for very low mass 
haloes, the accretion of subhaloes are perpendicular to 
63 . Regardness of the different approached, there are two 
possible additional reasons for the difference with our re¬ 
sults. Firstly, In Libeskind et al. (2014) they combine 
the signal for host halo mass normalized by the char¬ 
acteristic mass, M^,, of halo collapse. At 2 ; = 0, M* is 
about ^ 3 X 1 O^^M0. They found that for all haloes 
with m < O.IM* the accretion is along 63 . For a rough 
comparison, if we also combine the accretion signal for 
haloes with Mhost < O.IM* at z = 0, we also find that 
the alignment is along 63 . Secondly, the smooth length 
in Libeskind et al. is not constant for all haloes, but it 
is related to the virial radius of each halo. As we see in 
Fig. [5] using a small smooth length produces a stronger 
alignment between the halo major axis and the orienta¬ 
tion of density field around, thus it is expected that the 
alignment with the LSS is close to the alignment with 
the halo major axis, and being more universal. 

4. CONCLUSION AND DISCUSSION 

In this paper, we use N-body simulations to study the 
anisotropy of subhalo accretion with respect to the halo 
shape (major axis) and the large scale structure around 
(e 3 , the least compressed direction). We have obtained 
these main results: 

• The major axis of dark matter halo is well aligned 
with the least compressed direction, 63 , of the large 
scale structure around the halo, in agreement with 
the expectation from the Zel’dovich theory. 

• Subhaloes are accreted along the major axis of the 
host halo, and the alignment increases with host 
halo mass and redshift (stronger at high redshift). 
Most importantly, this alignment is universal (pos¬ 
itive) across all halo mass. 

• The accretion of subhalo along 63 is not universal. 
In massive host haloes, the accretion is aligned with 
63 , while in low-mass host haloes, the accretion of 
subhaloes is perpendicular to 63 . The transit mass 
for this flip is lower at high redshift. 

Our results show that there is a primordial anisotropy 
of subhaloes at accretion, and the degree of alignment is 
stronger than the observed alignment of satellites galax¬ 
ies at 2 : = 0. Thus, it is still difficult to conclude how 
much of the observed signal of galaxies is from the pri¬ 
mordial alignment or from the evolution effect inside dark 
matter haloes. In any case, the misalignment between 


central galaxies and host haloes is needed so as to ex¬ 
plain the observed degree of alignment between satellites 
and the shape of central galaxies (e.g., Dong et al. 2014). 

Our result about the accretion with the 63 direction 
has important implications on the correlation between 
the halo spin and the the LSS. Previous studies using 
N-body simulations (e.g., Zhang et al. 2015 and refer¬ 
ences therein) have reached convergent conclusion that 
the relation between halo spin and the LSS is not uni¬ 
versal. For massive haloes, their spin is parallel with 63 , 
while for low-mass haloes, their spin is perpendicular to 
63 . The characteristic mass for the flip of this correla¬ 
tion decreases with increasing redshift (e.g., Trowland et 
al. 2013). Some studies (Codis et al. 2012; Welker et 
al. 2014; Aragon-Calvo et al. 2014; Codis et al. 2015) 
pointed out that this is due to the disparity in the merger 
history of haloes with different mass. They speculated 
that for massive haloes, their subhaloes are accreted 
along the filament, and for low-mass haloes, their sub¬ 
haloes are accreted perpendicular to the filament. Un¬ 
der the assumption that the orbital angular momentum 
of accreted subhaloes is transferred into the host halo 
spin, it is naturally to expected a non-universality of the 
halo spin with filaments. However, such a spectulation 
for mass accretion disparity is not clearly shown in their 
analysis. Our results, for the first time, show that the 
accretion along 63 is really not universal, and the redshift 
dependence also agrees with the results from simulations. 

In deriving the above non-universality of the alignment 
between subhalo accretion angle with the orientation of 
the LSS ( 63 ), we use a constant smooth length for all 
haloes (although being evolving with reshift). Most pre¬ 
vious results (see Forero-Romero et al. 2014 for refer¬ 
ences) used a constant smooth length (ranging between 
0.5 — 2Mpc, well covered by our choice of smooth length), 
and most of them (if not all) reached agreement on the 
non-universality of the halo spin-LSS correlation. It im¬ 
plies that by using a constant smooth length for all halo 
mass, the non-universality on the accretion angle of sub¬ 
haloes must be in place so as to explain the simulation 
results. Our results do support this implication. Nev¬ 
ertheless, our choice of smooth length means that we 
are actually looking at the mass accretion pattern on a 
fixed scale at given redshift. As shown by some stud¬ 
ies (e.g., Aragon-Calvo et al. 2007a) using the multi¬ 
scale morphology filter to the density field, and they also 
found that the halo spin-LSS correlation has a mass de¬ 
pendence. This implies that the anisotropic accretion of 
subhalo could be both mass and redshift dependent, and 
it calls for a more complete study on subhalo accretion 
using different smooth length. 

In addition to the subhalo accretion anisotropy, to fur¬ 
ther understand in detail the correlation between halo 
spin and LSS, one needs to study the formation history 
of the halo spin in different environment of the LSS, in 
order to identify which kind of mergers contribute to the 
hnal correlation. Such an analysis should make use of the 
velocity and spin information from N-body simulations. 
Also, one has to understand why the correlation between 
subhalo accretion and 63 of LSS is different for high and 
low mass haloes. We will investigate these issues in a 
future paper. 
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